Additional carrier-mediated ferromagnetism in GdN 
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The mechanism behind ferromagnetic exchange interaction in GdN is not well understood. It has 
been argued that it can be due to fourth order cross process of d-f mixing and d-f exchange. An 
alternative explanation suggests an anti- ferromagnetic interaction between Gd d and N p induced 
moments on the rock salt structure which aligns the nearest neighbor Gd / moments ferromagnet- 
ically through the d-f exchange. In this paper we present results of Curie temperature in GdN 
as a function of carrier density calculated within our multiband modified RKKY- like exchange 
interaction. It includes realistic bandstructure of the 5d conduction band as an input for single 
particle energies. We analyze the possibility of carrier- mediated ferromagnetism in GdN and also 
demonstrate a simple phenomenological model which justifies the role of charge carriers. 

PACS numbers: 71.10.-w, 71.27. +a, 72.20.-i, 75.50.Pp 



I. INTRODUCTION 

Spintronics^ is a technology which utilizes the charge 
states of electrons as in a semiconductor as well as the 
quantum spin states as used in the data storage devices. 
The interest in semiconductor based spintronics has 
greatly accelerated the studies of magnetic semicon- 
ductors^. The aim is to find semiconductors exhibiting 
magnetism at or above room temperature. In order to 
achieve this a host semiconductor material was doped 
with transition metal ions (Mn) producing diluted mag- 
netic semiconductors (DMS) with Curie temperature, 
T c , as high as 185 Ki. Though there have been a lot 
of progress during last decades to accomplish this task 
but there is still lack of complete understanding due to 
complications in growth and measurement techniques 
and because of approximate theories. 

Among other class of materials are the strongly corre- 
lated rare earth (RE) compounds having incompletely 
filled / electron shells. In particular, Gadolinium Nitride 
(GdN) has been the most widely studied^ owing to its 
half- filled / shell with a magnetic moment of S=| which 
makes it more attractive for spintronics applications 
as compared to the transition metal doped materials. 
As far as its electronic properties are concerned, it 
was experimentally demonstrated to be a low carrier 
semi- metal 6 - and insulating-. There are also several 
recent reports 8 -^ of GdN having a degenerately doped 
semiconducting ground state based on the resistivity 
data measured at low temperatures. Theoretically it 
is predicted to have a semi- conducting 10,11 or a half- 
metallic character based on ah- initio calculation a 12 ' 13 . 

Initially there had also been a dispute regarding its 
magnetic properties with earlier observations describing 
GdN to be a metamagnelj 6 -^ (i.e., at low field an anti- 
ferromagnetic and at high field a ferromagnetic) material 
while other studies indicated it to be a ferromagnetic^—. 
However after such controversial discourse, it has been 
accepted that GdN is a ferromagnetio 8 ' 21 ' 22 material 
with experimental reported values of the T c in the range 
58 - 90 K-ii£r— . But there are as many puzzles on the 



mechanism of ferromagnetic exchange interaction in 
GdN as in case of its electronic ground state. 

This paper is organised as follows. In section |TT] 
we discuss the existing proposed mechanisms behind 
the magnetic ground state of GdN. In section IIIII we 
present results of our multiband modified RKKY- 
like exchange interaction which takes as an input for 
single particle energies the realistic bandstructure of 
the 5d conduction band. We examine the possibility 
of carrier- mediated ferromagnetism in GdN. In order 
to understand the origin of source carriers and their 
role in supporting ferromagnetism we consider a simple 
phenomenological model and study its validity. In sec- 
tionllVlwe summarize and conclude our obtained results. 



II. 



EXCHANGE INTERACTIONS 



GdN seem to behave like the ferromagnetic EuO 
since both have similar magnetic moments and values of 
T c . So one would expect the same mechanism behind 
exchange interaction in GdN as in EuO. But Kasuya 
and L i 24 : 25 explained the essential distinction which is 
briefly summarized below. 

There are three known mechanisms for the exchange 
interactions between magnetic atoms in rare- earth 
compounds. The first is due to the second order pertur- 
bation of the intra- atomic d-f exchange giving rise to 
RKKY 26 - like interaction. The second is due to third 
order perturbation theory of the d-f exchange and d-f 
mixing. The nearest neighbor exchange interaction Ji 
in Europium chalcogenides is due to this mechanism 27 . 
The third is due to the fourth order perturbation of 
the d-f or p-f mixings, where p are the anion states 2 ^. 
The first mechanism does not depend on the 4/ level. 
Whereas the second and third mechanisms become 
important when the 4/ level is near the Fermi edge as 
found in EuO. But in GdN, it is known 8 -^ that the 
4/ level lies much below the Fermi level. So the first 
could be one of the possible mechansims of exchange 
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interaction in GdN. 

Since both the compounds have different possible 
mechanisms then what could be the reason for having 
similar values of T c ?. Let us consider the nearest 
neighbor exchange interaction in Eu chalcogenides as 
briefly explained earlier and see how different it is in 
GdN. 

The exchange interaction, Ji is dominated by an 
indirect interaction arising from the virtual excitation 
of a 4/ (lying inside the semiconducting gap for EuO) 
to a 5d state, which then overlaps the neighboring Eu 
and leads to a /-/ interaction through the d-f exchange. 
This d-f exchange essentially measures the spin splitting 
of the d bands induced by their intra- atomic exchange 
interaction with the / state. One may visualize the effect 
as arising from the hopping of / electron to a neighboring 
site d orbital where it is subject to a spin exchange 
interaction, Jdf. In the language of perturbation theory 
it means that the d orbital gets mixed into / band in 
an amount where tdf is the hopping integral and 

(ea — ef) is the energy difference between the bottom of 
the d band and the localized / level. The contribution to 
the exchange interaction between the nearest neighbor / 
sites is inversely proportional to (ea — £f). 

One would expect a smaller Ji in GdN as the above 
mentioned energy difference is much larger in GdN 
with the 4/ level lying several eV below the conduction 
band edge as compared to EuO where it is inside the 
semiconducting band gap. Thus Ji in accordance with 
third order perturbation theory is an order of magnitude 
less in GdN than in EuO. So in order to have similar 
values of T c in both the compounds there should be 
another dominating indirect exchange mechanism in 
GdN. 

Kasuya and L i 24 ' 25 developed the fourth order per- 
turbation theory which considers the cross process 
between the d-f mixing and d-f exchange interaction. 
The resulting effective spin- spin exchange interation 
depends on inverse of the energy gap i.e., difference 
between the bottom of the d band and the top of the 
p band, (ea — e p ). This exchange interaction energy is 
evaluated to be large in GdN because the theoretically 
ascribed energy gap in GdNi£ is small as compared to 
the gap in EuO^. But since the experimental nature 
of the electronic ground state in GdN is not yet clear 
and the value of the energy gap has not been reported 
this mechanism remains to be verified. And moreover as 
mentioned earlier that the 4/ level in GdN is several eVs 
below the Fermi level so the mechanisms due to fourth 
and third order seems less plausible. 

Recently, Mitra and Lambrecht^ presented an alter- 
native way to explain the ferromagnetic ground state 
structure in GdN. There is as an anti- ferromagnetic 
ordering on a rocksalt lattice between N p and Gd 
d magnetic moments. And due to the d-f exchange 
coupling, the nearest neighbor Gd atoms interact ferro- 
magnetically with each other. According to their picture 
even the next nearest neighbors are ferromagnetically 



aligned. But they obtained a T c of 10 K within a mean- 
field calculation which is much lower than observed 
experimentally. 

So is there an additional indirect exchange mechanism 
in GdN ? In recent experiments^, GdN films were found 
to be semiconducting doped to degeneracy with the most 
likely source of charge carriers (electrons) as nitrogen 
vacancies. Although they found the same value of T c 
as before and also the Hall effect measurements showed 
the presence of charge carriers but their study with 
further doping of carriers suggested that exchange is 
not mediated by free carriers. And this view support 
the earlier theoretical explanation 3 -^ where the authors 
abandoned the existence of RKKY- like exchange within 
rigid band model. 



III. MODELS, RESULTS & DISCUSSION 

Here we consider the carrier- mediated ferromagnetism 
in GdN as a possible additional exchange interaction. 
But our theory is unlike the one studied earlie r 31 ' 32 . We 
are interested in determining the magnetic properties of 
multiband Kondo lattice model which is described by the 
following Hamiltonian, 

H = Hkin + Hi n t (1) 

where 

H ki „ = £ Tf cLc j/3(T (2) 

ijapcr 

is the kinetic energy of the system and 

Hint = — ^( z <rS^C i f Q0 .Ci Q0 - + Sf cJ Q _ CT Ci QO -) (3) 

is the intra- atomic exchange interaction term with an as- 
sumption that itinerant electron in each band is coupled 
to the localized moment by the same coupling strength, 
Jdf. The Greek letters (a,/3) depict the band indices. In 
GdN, these are the five d conduction bands. The latin 
letters (i,j) symbolize the crystal lattice sites and spin is 
denoted by a{— t,i)- 

The total Hamiltonian, Eq. (fT]), can be solved using 
Green function method. Since our interest lies in mag- 
netic properties of multiband Kondo lattice model we 
have to consider both the sub- systems (localized as well 
as itinerant) within a self consistent scheme. In Ref^ we 
have presented our modified RKKY theory which treats 
both the sub- systems equally. The main idea of the mod- 
ified RKKY theory is to transform the above Kondo- like 
exchange Hamiltonian of the conduction electrons into 
an effective Heisenberg- like spin- spin exchange Hamil- 
tonian of the / spins by averaging H; n t in the subspace 
of the conduction electrons. In order to avoid repetition 
of any kind we refer to the reader Ref £^ for the complete 
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analysis. But in the following we would like to highlight 
the numerical details on calculating the exchange inte- 
grals and thereby the Curie temperature. 

The sensitivity of the RKKY- like mechanism to car- 
rier concentration is well known and that the RKKY os- 
cillations are strongly dependent on the value of carrier 
concentration. The distinct method&-&21 ~ 23 i 33 of sample 
preparation lead to different carrier concentrations. This 
could not only be reason why GdN was noted to be an 
anti- ferromagnet at low fields with Neel temperature of 
40 K— iM but also why the experimental determination of 
electronic ground state (semi- metallic or semiconduct- 
ing) is so uncertain. 

We take into consideration the semi- conducting na- 
ture as predicted theoretically^ and obtained experi- 
mentally^. In the latter case, the carrier concentration 
(doping) is usually assigned to defects like nitrogen va- 
cancies^ or structural defects (grain boundaries between 
the nanocrystallites)^ 3 ". With realistic values of input pa- 
rameters^ like strength of d-f exchange coupling and the 
single particle energies of 5d conduction band obtained 
using TB-LMTO-ASA 36 we evaluate the following effec- 
tive exchange integrals, 



J cff (q) 
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(5) 

_ — (7 

In above equations, G (E) and G (E) are the sin- 
gle particle non- interacting (undressed) and interacting 
(dressed) Green function matrices respectively and f_ (E) 
is the Fermi functional. The subscript Y denotes the s th 
neighboring shell of radius R s spanning As number of 
neighbors to the central atom as shown in the inset of 
Fig. [T] for the case of planar geometry and only for first 
six shells. But this notation can be generalized to any 
lattice and finite number of shells until convergence for 
the exchange integral is reached. 

In Fig. [1] we also present the dependence of effective 
exchange integral on the distance for a few low values 
of carrier concentration. In our numerical calculations 
we have considered upto 40 shells in 3D geometry of fee 
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FIG. 1: (Color online) The effective exchange integral as a 
function of distance for various values of band occupation. 
The first six shells of nearest neighbors to a central point are 
shown in a planar geometry. 



structure for the Gd atoms. As seen in Fig. [TJ the os- 
cillations damp out as a function of distance giving a 
typical characteristic long range RKKY- like behavior. 
The atypical part is that the strength of next nearest 
neighbor interaction is stronger than the nearest neigh- 
bor. But it keeps on decreasing rapidly as we increase 
the carrier concentration and eventually becomes weaker 
than the nearest neighbor for n=0.1. 

As the Curie temperature is dependent on effective ex- 
change integrals in the following form 3 !, 



T r = 



2S(S+1) 
3k B 



1 



J cn (0)-J ctt (q)7 T 



(6) 



where J cff (0) = J cff (q = 0), we calculate T c within a self- 
consistent scheme for the values of carrier concentrations 
as considered earlier. 

Fig.[2]exhibits the dependence of T c on carrier concen- 
tration for GdN with the experimentally reported range 
depicted within the horizontal lines. We obtain our high- 
est value of T c =60.144 K for a carrier concentration, 
n=0.01 (~ 8 x 10 19 /cm 3 for lattice constant of 5 A). 
And for concentration of n=0.1 (~ 8 x I0 20 /cm 3 ) the 
T c drops down to zero. Such high values of carrier con- 
centration have been observed experimentally 9 -. 

In our theory for each value of carrier concentration we 
determined the Fermi edge self- consistently. We didn't 
regard any impurity level or band in our calculations. It 
would be interesting to model the source of charge carri- 
ers in a realistic way. And study the contribution of such 
defects or impurities in stabilizing the carrier- mediated 
ferromagnetism in GdN. We keep this analysis for future 
work. But in order to confirm the role of charge carriers 
and explain one of the anomaly in the low temperature 
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FIG. 2: (Color online) The dependence of T c on band occupa- 
tion as calculated using modified RKKY theory for Jdf=0.353 
eV and S=3.5. The reported experimental range is within 58 
- 90 K as shown by two horizontal lines. The dotted black 
line is a guideline for the eye. 



behavior of resistivity^, we examine a very simple phe- 
nomenological model as shown in Figure [3] 

We consider the spin splitting of the parabolic model 
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FIG. 3: (Color online) A schematic depicting the spin split- 
ting of the density of states below and at the T c . A movable 
impurity level (dashed straight line, E=Eo) is also shown with 
a thick arrow pointing its direction. 

density of states below T c within the mean field picture 
having energy difference of the order ~ J^f S where S is the 
4/ local moment. An impurity level which can be moved 
in or out of the spin-up bands is also shown. Although re- 
cent experiments 9 support the presence of charge carriers 
(electrons) due to impurities but no evidence is reported 
on whether do these carriers form an impurity level or 
degenerate impurity band. We take into account an im- 



purity level instead of a degenerate band. It is known 
that the impurity concentration in case of intrinsic semi- 

— AE 

conductor at room temperature is proportional to e k B T . 

AE 

So we consider the resistivity to be proportional to e k B T 
and is given as, 



p = p(0) e k B T 
where the activation energy, AE, is given by 

AE=(E g -Eo) + ^ (SZ) 



(7) 



1-V (8) 



The reason to assume such a form of activation energy 
can be understood as follows. Our goal is to model 
a ferromagnetic semiconductor with a finite gap E g 
3> keT and described within Kondo lattice model. It 
governs the temperature dependence via the magneti- 
zation ((S z )) which we consider within the molecular 
field theory. Below the ferromagnetic T c the resistivity 
depends on the scattering of charge carriers due to their 
interaction with localized moments. Above T c , the 
resistivity follows the normal thermally activated energy 
behavior and falls off exponentially with further increase 
in the temperature. 

The first term in the bracket of Eq. ((SJ represents 
the loss in binding (trap) energy of the electrons from 
the impurity level and the second term is the exchange 
contribution that reaches its maximum value at or above 
T c . Using the parameters from the literature^ ' 10 ' 35 i.e. 
J df =0.353 eV, S=| , E g =0.30 eV, T c *p = 68 K and 
p(0)=0.4 Ocm, we try to evaluate Eq. © for different 
values of Eo. We move the impurity level from the center 
of the gap into the lower edge of the spin- up density of 
states of the conduction band. The result is shown in 
Figure |U 

Our simple model demonstrates an important effect 
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FIG. 4: (Color online) Temperature dependent resistivity ex- 
hibiting metal- to- insulator transition upon the variation of 
the impurity level, Eo, out and in the lower edge of conduction 
band. 

like insulator- to- metal transition at low temperature by 
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moving the level Eo as shown in Figure [3] Similar effect 
induced by applying strain was theoretically reported 
based on first principles calculations 3 - 7 .. Thus it clearly 
reveals that the experimentally observed resistivity^ is 
an artifact of a degenerate impurity level lying close to 
the lower edge of d conduction band. And carriers in 
such level play an essential (additional) part in stiffening 
the ferromagnetism in GdN. 



IV. SUMMARY & CONCLUSION 

In conclusion the basic exchange mechanism in GdN 
is not clear. There is a theory based on fourth order 
perturbation^ which seem less probable since it requires 
the 4/ level to be near the Fermi edge while experi- 
mentally it is known to be several eVs below. There is 
another 30 interpretation but it results in obtaining much 
lower T c as compared to the experimentally reported 
values. 

We consider our multiband modified RKKY theory 
where an effective spin Hamiltonian is obtained by 



integrating out charge degress of freedom from the 
multiband Kondo lattice model^i. We take the realistic 
bandstructure of 5 c? conduction bands as an input 
for the single particle energies and the d-f exchange 
coupling 35 to calculate the dependence of T c on carrier 
concentration. The results are in close proximity to the 
experimental findings. 

In order to trace the source of charge carriers which 
are eventually responsible for an additional carrier- 
mediated ferromagnetism in GdN we consider a simple 
phenomenological model. It not only explains exper- 
imentally observed anomaly in the low temperature 
behavior of resistivity^ but also exhibits insulator- to- 
metal transition in accordance with theory 3 — Our 
results indicate that if pure (stoichiometric) GdN is 
prepared it will have a low T c as predicted earlier—. 
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